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X-ray flares from Sgr A*: star-disk interactions? 
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Abstract. Sgr A*, the putative black hole in our Galactic Center (GC), is extraordinary dim in all frequencies. 
Apparently the black hole is unable to accrete at the Bondi accretion rate for some reason. Another mystery of 
Sgr A* is the recently discovered large magnitude and short duration X-ray flares (Baganoff et al. 2001). Here we 
point out that X-ray flares should be expected from star passages through an inactive (i.e. formerly accreting) 
disk. There are thousands of stars in Sgr A* stellar cluster, and each star will pass through the disk twice per 
orbit. A shock hot enough to emit X-rays is formed around the star when the latter is inside the disk. We develop a 
simple yet somewhat detailed model for the X-ray emission from the shock. The duration of the flares, their X-ray 
spectra, frequency of events, and weakness of emission in the radio and near infra-red appear to be consistent 
with the available observations of X-ray flares. We therefore suggest that at least a fraction of the observed flares 
is due to the star-disk passages. Such star-disk flares are also of interest in the nuclei of nearby inactive galaxies, 
especially in connection with perspective Constellation-X and XEUS X-ray missions. 
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1. Introduction 

1.1. Sgr A* : underfed X-ray flaring black hole 

With the current high quality data (e.g., Schodel et al. 
2002, Genzel et al. 2003; Ghez et al. 2003a,b) it seems very 
safe to say that the center of our Galaxy hosts a super- 
massive black hole, with Mbh — 3 x 10^ Mq. The black 
hole is associated with Sgr A* source, the compact ra- 
dio core (Rcid ct al. 1999), discovered by Balick & Brown 
(1974). Thanks to Chandra observations (Baganoff et al 
2003a), we know that there is enough of hot gas at dis- 
tances of the order ~ 10^^ cm (or ~ 1") to make the black 
hole glow at around ~ few xlO'*° erg/sec if accretion pro- 
ceeds via the standard radiatively efficient disk (Shakura 
& Sunyaev 1973). However the observed bolometric lumi- 
nosity of Sgr A* is 4 orders of magnitude smaller than 
expected (e.g. Zhao, Bower & Goss 2001; see also review 
by Melia & Falcke 2001). The luminosity in X-rays is 
even smaller and is ^ 10"^ of the nominally expected 
numbers (Baganoff et al. 2003a). Understanding of this 
discrepancy is the key to accretion physics at low accre- 
tion rates which encompasses a very large scope of objects 
from galactic centers to stellar mass compact objects (e.g., 
Narayan 2002). 
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The currently leading explanation of the low luminos- 
ity from Sgr A* is the Non-Radiative Accretion Flows 
(NRAFs), which were discussed in great detail first by 
Narayan & Yi (1994). The corner stone assumption of this 
model - that ions are much hotter than the electrons in 
the inner region of the flow - has not yet been indepen- 
dently verified. Initially it was believed that these flows 
carry ("advect") most of their energy into the black hole, 
radiating very little on the way. Several years ago it be- 
came clear that these hot flows also drive powerful thermal 
winds (Blandford & Begelman 1999), so the accretion rate 
onto the black hole is significantly reduced. This theoret- 
ical picture seems to be confirmed in its gross features by 
the observations (e.g., Narayan, Yi & Mahadevan 1995; 
Bower et al. 2003; Quataert 2003), although alternative 
interpretations of Sgr A* also exist (Melia 1994; Falcke 
& Markoff 2000; Coker & MeUa 2000; Yuan, Markoff & 
Falcke 2002). 

The recently discovered X-ray flares in Sgr A* 
(Baganoff et al. 2001) brought further theoretical chal- 
lenges. The X-ray flares happen at a rate of about one per 
day (Goldwurm et al. 2003a; Baganoff et al. 2003b). The 
X-ray luminosity during the flares exceeds the quiescent 
one by up to a factor of 100 (although most of flares are 
weaker than that). Flare durations range from ^ 1 to 10 
kiloseconds. The X-ray spectrum of flares is harder than 



2 



Nayakshin et al.: X-ray Flares and inactive disk in Sgr A' 



the quiescent one. The follow-up multi-wavelength cam- 
paign of Sgr A* has not detected significant variability in 
any other band except for X-rays (Baganoff et al. 2003b; 
see also Hornstein et al. 2002). This latter fact is partic- 
ularly troubling since in all the current models the X-ray 
flares occur at the inner part of the flow that emits Sgr A* 
radio emission. Thus the challenge is to understand why 
such a large variation in X-rays does not lead to flares in 
other frequencies. 

1.2. Cold disks in LLAGN and Sgr A* (?) 

It is believed that accretion flow physics at low accretion 
rates should be largely self-similar (e.g. Narayan 2002) and 
the difference is just in scales (e.g Mbh and accretion rate 
M). Therefore Sgr A* is expected to be physically simi- 
lar to the Low Luminosity AGN (LLAGN; e.g. Ho 1999). 
Most of LLAGN appear to posses thin cold disks. The 
presence of these disks is deduced convincingly through (i) 
water maser emission (Miyoshi et al. 1995); (ii) infrared 
"bumps" in the spectral energy distribution of LLAGN 
(Ho 1999) that are thought (Quataert et al. 1999) to be the 
thermal disk emission (with disk temperatures T ^ 10^ 
K); and (iii) the prevalence of the double-peaked broad 
emission lines (e.g. Ho 2003). The outer radii of these disks 
are poorly constrained, whereas the inner disk radii can 
be constrained from both the SED and the emission line 
profiles, and seem to be in the range ^ 100 — lOOOi?^, 
where Rg is the gravitational radius. 

It is thus tempting to suggest that such a disk may 
also exist in our GC (see t l2.1|l . There appears to be sev- 
eral promising ways to test the disk hypothesis. The sim- 
plest of all is through eclipses of close stars (Nayakshin & 
Sunyaev 2003). In particular, the orbit of star named S2 
is very well determined (Schodel et al. 2002, Ghez et al. 
2003b). Another important effect is the reprocessing of the 
optical-UV radiation of the star in the disk into the near 
infrared (NIR) band, which can yield a strong increase in 
the perceived NIR flux of the star. The current S2 data 
alone practically rule out an optically thick disk (Cuadra, 
Nayakshin & Sunyaev 2003). However a disk optically thin 
in NIR is permitted. 

1.3. X-ray flares in Sgr A* : star disk interactions? 

In this paper we point out an additional effect as a test of 
the inactive disk hypothesis. Recent results of Gcnzel et 
al. (2003) show that there are as many as ~ lO'' stars in 
the inner arcsecond lO^Rg) of Sgr A* . Several of these 
stars cross the inactive disk every day. While inside the 
disk, the stars drive a very strong and hot shock wave in 
the disk material. The shock will emit X-rays (since the 
shock is short lived, this emission is actually a flare). In 
this paper we shall build a simple yet somewhat detailed 
and self-contained model for these events. 

Let us point out that such star-disk flares have a po- 
tential to naturally explain why there are no counterparts 



of X-ray flares in the radio and near infrared frequencies. 
A somewhat unexpected feature of this model is that the 
flares are the result of a process that has nothing to do 
with the hot gas that is responsible for Sgr A* quiescent 
emission in all frequencies. This is so because stars are 
"small" (~ 10^^ cm) compared with the size of the hot 
flow (~ 10^^ cm), and they have a negligible direct effect 
on the flow of hot gas. Flares are emitted when the stars 
strike the cold disk that is much denser than the hot flow. 
This happens at radii of order lO'^ — lO^Rg, and thus the 
flares do not interfere with the flow at i? ;S lOg where the 
jet originates. In addition, the X-ray flares are emitted 
by a thermal moderately hot gas, and we will later see 
that the radio and near infrared flares from the star-disk 
passages are weak. 

Finally, the problem of X-ray flares from star-disk pas- 
sages merits a careful study even if there is no disk in our 
own Galaxy. Similar X-ray flares may be detectable in 
other LLAGN, especially the close ones such as M31 or 
M32 (Ho, Terashima & Ulvestad 2003). 

The plan of the paper is as following. We flrst attempt 
to determine in ^ the kind of a disk that may reside in 
Sgr A* from arguments that are independent from the ob- 
served X-ray flares. We flnd that the disk must be inactive, 
i.e. not accreting, with temperature ^ few xlO^ K, with 
gas density in the range of 10^" to 10^^ cm""^. In H2.2I we 
will see that a disk optically thick in near infrared is ruled 
out but an optically thin disk is allowed. In l^l'^c calcu- 
late the rate or the star-disk passages, and also show that 
the disk is very likely to have an inner hole. In 21 we show 
schematics of a star-disk passage and mention the connec- 
tion of our work to the existing literature on such events. 
X-ray luminosity and duration of flares is calculated in 
SjSl X-ray spectra are considered in SjHJ whereas the near 
infrared spectra are addressed in ^0 In ^properties of 
flares produced by a low mass star and a high mass star 
are discussed, and in ^l^a comparison between the model 
and observations is made. A list of shortcomings of our 
work is given in §refsec: short. Finally, discussion of the 
results is presented in ijlll 



2. Constraints on the disk properties 

2.1. Arguments for and against a disk in Sgr A* 

If the putative accretion disk were of the standard type 
(Shakura & Sunyaev 1973), then it would have been seen 
by now even if accretion rate is as small as 10~^ Mo/year 
(Narayan 2002). But at very low disk accretion rates hy- 
drogen in the disk becomes neutral and the accretion rate 
plunges to very low values (e.g. see review by Cannizzo 
1998; see also equation El below). This is why these disks 
are inactive, i.e. not accreting. In addition, these disks 
seem to miss their inner regions for a variety of reasons 
(e.g., Quataert et al. 1999 and i j3.2() and therefore they 
can escape the near infra-red (NIR) detection limits that 
rule out standard steady-state disks. 
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However, Falcke & Melia (1997; FM97 hereafter) 
probed the existence of the disk through its interaction 
with the winds from hot stars. This wind is beheved to 
be the source of the hot gas that should be accreting onto 
the black hole. FM97 assumed that the wind runs into 
and settles down onto the disk, releasing the mechani- 
cal energy into the black-body radiation on impact. They 
found that unless the wind has a large specific angular 
momentum (so that it cannot penetrate to small radii), 
the luminosity expected in the NIR violates the available 
limits on Sgr A* NIR emission. Since there are ^ 10 wind 
sources, no large angular momentum in the wind was ex- 
pected. FM97 therefore stated that existence of the disk in 
Sgr A* is highly unlikely (see also Coker, Melia & Falcke 
1999). 

New data show that eleven out of twelve wind- 
producing stars rotate in the same (counter-Galaxy) di- 
rection (Genzel 2000). These and other hot luminous stars 
in the region share the same puzzling property: they are 
quite young (Ghez et al. 2003b; Genzel et al. 2003) and it 
is unclear how these stars made it to the innermost arc- 
second of the Galaxy. One of the possibilities is that these 
stars were created from a single large molecular cloud that 
fell into the inner region of the Galaxy with a small impact 
parameter. This cloud could have created both the stars 
and an accretion disk whose (tiny not accreted) remnant 
may still be there. ^ Note also that matter in the accre- 
tion disk flows inward and outward from the characteristic 
(circularization) radius (Kolykhalov & Sunyaev 1980). It 
is thus not possible to accrete all of the disk in the active 
phase without leaving a remnant. 

Nayakshin (2003) recently showed that, in the setting 
similar to FM97, the hot flow rapidly looses its heat to 
the inactive disk. The latter re-radiates this energy in the 
infrared. As the hot gas looses heat, it also looses vertical 
pressure support, and as a result condenses onto the inac- 
tive disk. The accretion flow in the radial direction is es- 
sentially terminated. Thus the hot gas is accreting not on 
the black hole but rather on the inactive disk. The tightest 
of the existing NIR limits is not violated by the model if 
wind circularization radius is i?c ^ 3 x 10'*i?g. Since the 
wind sources are a few arcsecond few x lO^i?^) away 
from Sgr A* (e.g. Genzel et al. 2003), it appears that this 
"large" value of the circularization radius is quite reason- 
able. The luminosity emitted by Sgr A* is very small in 
this model because the gravitational potential at the con- 
densation point is tiny compared with that of the last 
stable orbit around the black hole. We refer to this model 
as the "Frozen Accretion" (FA) model of Sgr A*, and use 
some of the constraints from the model to limit the pos- 
sible size and density in the disk. 



2.2. Disk temperature (and NIR optical depth) 

First consider a disk optically thick in all frequencies. Let 
Rd be the outermost radius of the disk. Since the disk is 
inactive, the effective temperature of such a disk does not 
behave as the usual oc For example, appears 

to stay nearly independent of radius for quiescent disks 
of eclipsing dwarf novae (e.g. Figs. 5 & 6 in Bobinger et 
al. 1997). If Td ~ const, then the disk total energy out- 
put is dominated by the emission from outermost radii. 
Accordingly the disk luminosity is Ld ^ 2iTR^aBT^. In 
order to satisfy observational constraints on the total lu- 
minosity of Sgr A* in quiescence, Lq, this temperature 
must be smaller than 



in 



2t:R^ con 1(7 b 



1/4 



= 77K 



.1/4 
^36 



■1/4 1/2 



(1) 



where L36 ~ Lq/10^^ erg/sec, r4 = Rd/ 10 Rg, Rg = 
2GM/c^ = 9 X 10^^ cm is the gravitational radius for 
the black hole mass of Mbh = 3 x 10^ Mq. We shall see 
below that the disk is unlikely to be very optically thick 
so the disk midplane temperature should not exceed the 
effective temperature by a large number. 

We should also require the disk to be hot enough to 
re-radiate the visible and UV radiation of the Sgr A* star 
cluster itself. We use the latest results on the star number 
counts in Sgr A* cluster from Genzel et al. 2003, together 
with the fact that the total luminosity of stars in the cen- 
tral parsec of the Galaxy is about 50 million Solar lumi- 
nosities (^0). We first calculate the light to mass ratio, 
and then use equation El to find that the integrated light 
of the star cluster within radius R is 



L^{R) = 2.1 X 10" Lq rl''^ erg/sec 



(2) 



The flux incident on the disk is then F^,{R) — L^,/4ttR'^, 
and the resulting effective temperature due to stellar ra- 
diation heating is 



n = 190 K r, 



-1/8 



(3) 



which is a factor of 3 higher than the value in equation ^ 
This latest result shows that an optically thick disk in fact 
would over-produce the infrared luminosity of Sgr A* by 
about a factor of a hundred (e.g. see Sgr A* spectrum in 
Falcke & Melia 2001 and in Narayan 2002). Therefore the 
disk must be optically thin in the infrared frequencies or 
highly inclined with respect to our line of sight. An esti- 
mate shows that for not too small disk inclination angles 
the infrared optical depth of the disk should be no larger 
than about 0.01 at frequency ~ 10^^ Hz. In addition, close 
passages of very bright stars such as S2 (L ^ 10^ Lq; Ghez 
et al. 2003b) may actually increase the value of to about 
300 K (Cuadra et al. 2003). 



^ We have recently learned that, using data of R. Genzel's 
group. Levin & Beloborodov 2003 found that 10 out of 13 wind- 
producing stars may lie "exactly" in a single plane, which fur- 
ther confirms our suggestions. 



2.3. Disk size 

The quiescent spectrum of Sgr A* is in rough accord with 
the data if the main part of the hot flow condenses onto the 
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cold disk at radius i?c ^ 3 x 10^ Rg (Nayakshin 2003). If 
condensation occurs at much smaller radii, then the flow 
gains too much gravitational energy (while falling from 
the Bondi radius, Rb ^ IQ^Rg to Rc) and the luminosity 
of the disk would be too high. The value of Rc could be 
however reduced if the exact accretion rate in the con- 
densing flow is lower than the Bondi estimate (similar to 
the ADAF case; see §5.2 in Narayan 2002). Therefore we 
estimate the minimum disk outer radius as 

Rd-^l-lOx 10*i?g (4) 

There are no strong constraints based on the FA model on 
the inner disk radius, Ri, as long as it is at least a factor 
of few smaller than condensation radius Rc- 



photo-ionized layer of gas, and compared it with that of 
the S2 star in the K-band. This flux is about 10^'^ of the 
star's continuum vFnu at that frequency. For comparison, 
this is about 30 times weaker than the absorption in the Br 
7 line from the star itself (see Fig. 1 in Ghez et al. 2003b). 
In addition, the emission line from the photo-ionized layer 
should be 5-10 times broader than the line from the star 
due to much faster disk rotation, making it challenging to 
observe. 

Finally, free-free absorption in the photo-ionized layer 
is important for frequencies smaller than ^ lO^^rj^,^^ Hz. 
For these frequencies, the corresponding thermal radio 
emission of the layer is at the level of about 1 percent 
of the Sgr A* jet emission at the respective frequencies, 
rendering it negligible. 



2.4. Constraints from close bright star passages 2.5. l\/1inimum disk density 



The second string of constraints on the disk size and in- 
clination angle, i, comes from the fact that close stars 
could be eclipsed by the disk if it is optically thick in 
near infrared. Here we summarize our preliminary results 
(Cuadra et al. 2003; see also Nayakshin & Sunyaev 2003) 
based on the measured S2 star's positions. The parame- 
ters of the orbit are now very well known (Schodel et al. 
2002; Ghez et al. 2003b). Cuadra et al. (2003) find that if 
the disk is optically thick and has an inner hole less than 
^ 0.01" ~ 10'^ Rg, then the outer disk radius may not 
much exceed 10^ i?,,, barely making the constraints based 
on the FA model. In addition, for an inner hole of ~ 0.03 
arcsecond, there is a large range in the disk orientation 
parameter space that no eclipses of S2 are predicted for 
any arbitrarily large value of the outer disk radius. 

However reprocessing of the optical-UV radiation re- 
ceived by the disk from the star into the NIR band could 
significantly increase the apparent star brightness if the 
star is close enough to the disk. Using the deduced S2 
orbit, we studied this effect for all possible disk orienta- 
tions, assuming an optically thick disk with no inner hole 
and infinitely large outer radius. We find (Cuadra et al. 
2003) that an increase in the 2.2/im flux of S2 should have 
been detected at least in one of the S2 measured positions, 
whereas observations do not confirm that. Therefore the 
disk, if any, must be optically thin in the NIR. We are 
currently modeling similar reprocessing for an optically 
thin disk. Unfortunately these constraints are not easily 
converted into constraints on the total column depth of 
the disk, E (see t|2.8|l . because the properties of the dust 
grains in the disk are very uncertain (see 

We have also checked some other potentially observ- 
able indicators of the disk presence. Due to star's UV ir- 
radiation, a layer of completely ionized hydrogen forms 
on the top of the disk. Assuming temperature of order 
10"' K for the layer, we found that its column depth is 
Nh ~ lO^^niir^^ cm^^, where ris is distance between 
the star and the disk in units of lO"'^ cm, and rin is the 
hydrogen nuclei density in units of 10^^ cm~'^. We then 
estimated the Brackett 7 (A = 2.16/im) line flux from this 



Pressure equilibrium argument. In the frozen accretion 
model of Nayakshin (2003), the inactive disk and the hot 
accretion flow are in direct physical contact. Thus the disk 
midplane pressure should be at least as large as the pres- 
sure in the hot flow. Assuming the disk consists mainly 
of molecular hydrogen, Pd ~ fc7id2d/2, where rid is the 
midplane density of hydrogen nuclei in cm^'^. The den- 
sity in the hot accretion flow can be estimated through 
M ^ A-KR^VRmprihot, where vr ~ a^/ GMbh/R and 
a — O.lai is the viscosity parameter (Shakura & Sunyaev 
1973). This yields 



Tid > 1.5 X 10^ cm 



M'T' 
aiT2 



3 X 10'*i?„^ 



Rc 



(5) 



where M' = M/3 x 10"^ Mq year~^ is the accretion rate 
in the hot flow, T' = rhot/3 x 10^ K is dimensionless tem- 
perature of the hot flow, Rc is the wind circularization 
radius, and T2 — Td/100 K is the disk unperturbed mid- 
plane temperature. 

Total accreted mass argument. If condensation of the hot 
flow has been going on for a time ta, then the mass ac- 
creted by the disk from the wind is AM ^ Mta. If the 
accreted mass is deposited over area ^ Tri?^, then the ac- 
creted surface density is 

AE, ~ 3 g/cm^ Af'U ( ^^1^^) ' , (6) 



where tae = years. This translates into the disk 

density increase of 

And ~ 4 X 10^° cm-^ M'ta^T^'^' ( iiii^'^^' 

V Rr 

(7) 

This estimate yields a larger value for the minimum disk 
density than does equation ||SJ), so we conclude that the 
minimum disk density is in the range of ^ few x lO""' 
hydrogen nuclei in cm'^. Note however that the above ar- 
guments are applicable to the disk radii with R Rc, and 
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unfortunately they cannot be robustly extended to radii 
much smaller than that. 



2.6. Disk thickness 

The disk half thickness, H, is given by the hydrostatic 
balance (we neglect self-gravity; see below): 



We can also calculate the accretion rate through the 
disk as Mj, ^ 2i:RY.aCs{H / R)"^ , which yields 



10 



-10 



year 



7/2^ 

aniir/ 1. 



2 • 



(13) 



H = 



GMm„ 



= 3.9 X 10^2 T^'^rV^ cm 



(8) 



Note that (i) this accretion rate is extremely small; (ii) it 
is not a constant with radius since the disk is not a steady 
state one and of course in a real disk density rin will 
vary with radius. Nevertheless this equation clearly shows 
that even if viscosity parameter a is reasonably large, i.e. 
a ^ 0.1, the disk is extremely dim. 



where we assumed that the gas is dominated by molec- 
ular hydrogen so /i ~ 2mp. Note that the star's ra- 
dius (~ 10'^ cm) is much smaller than the disk height 2 9 Summary of disk properties 



scale for R > lO'^i?^. The disk is very flat since H/R 
0.96 X IQ-^T^^^ rl^^ . 



2. 7. Mass of the disk and maximum density 

The mass of the accretion disk with a constant (with ra- 
dius) midplane density and temperature, is 



Md~ 0.1 Menu 



/ Rd 



\ lO^Rn 



7/2 



1/2 



(9) 



where rin = nd/10 cm~ is the midplane density of hy- 
drogen nuclei. Md is much smaller than the black hole 
mass so that the disk is not globally self-gravitating 
(i.e., in the radial direction; Shlosman & Begelman 1989) 
and will not gravitationally influence stellar dynamics in 
Sgr A*. The disk is locally (vertically) self-gravitating 
when Q = Md/{H/R)MBH ^ 1 (Goldreich & Lynden- 
Bell 1965; Kolykhalov & Sunyaev 1980; Gammie 2001). 
We have 



= 5.0 X 10"^ r|nii , 



(10) 



iH/R)MBH 

i.e. the disk is not vertically self-gravitating as well for 
the nominal values of rid and Rd in this equation. We 
can turn the question around to derive the maximum disk 
density as the density at which the disk would become self- 
gravitating. In the latter case it would presumably break 
up into self-gravitating rings or regions and the star-disk 
impacts would become too rare. Hence, 

£ lO'^rf (11) 

For r4 = 3 equation (|11|) yields maximum density 
^ few xlO^'^ cm~3. Note that equation 1111 is applicable 
for any radius in contrast to equations |31 and 



Using quiescent Sgr A* spectral constraints, we found that 
the disk temperature should be Td ~ 100 K. The disk den- 
sity should be small enough to preclude it from being self- 
gravitating and large enough to accommodate the mass 
deposition from the hot flow settling onto the disk. The 
respective maximum and minimum values of the disk den- 
sity are (according to equations 171 & [TT|l are ^ few x lO"'^^ 
and few x 10^" cm^^ at radius i? = 3 x 10"^ Rg. If wind 
circularization radius is Rc = W^Rg then these values in- 
crease to ~ 10"'^^ and ~ 2 x 10^^ cm"^, respectively. With 
this range in disk density and also possibly large dynam- 
ical range between the inner and outer radii the disk can 
be either optically thin or thick in both X-rays and NIR. 
If the disk is optically thick in NIR, then most likely there 
must be a relatively large inner hole {Ri ^ 3 x lO^i?^) in 
the disk or else the disk would have eclipsed (Cuadra et 
al. 2003; Nayakshin & Sunyaev 2003) the S2 star (Schodel 
et al. 2002; Ghez et al. 2003b). However most likely the 
disk is optically thin in the infrared frequencies (see ^2.2\\ . 

3. The role of the star cluster 

3.1. The rate of star-disk crossings 

The distribution of stars in the central sub-parsec region of 
the Galaxy has been a subject of intense research for many 
years (see review by Genzel 2000). Existence of a stellar 
cusp with density distribution p^,{R) oc R^^-^^ has been 
predicted by Bahcall & Wolf (1976). Alexander (1999) 
found that a power-law stellar density distribution with 
exponent p between 1.5 to 1.75 is statistically slightly fa- 
vored over a constant density one. With unprecedentedly 
high statistic, Genzel et al. (2003) determined the stellar 
density distribution to follow the law 



2.8. Column depth and accretion rate 

The disk surface density is S ~ 2Hndmp, and is equal to 



S = 1.3 niir^^^ T2^^g cm 



(12) 



It is therefore clear that with the minimum and maximum 
disk density values just found, the disk can be either op- 
tically thin or optically thick to X-rays, and we should 
explore both cases. 



p*iR) 



p,o{R/Rbyp, R<Rb 
p^,a{R/ Rb)^^ , Rb < R 



(14) 



where i?b = 10 arcsecond,p = 1.4±0.1 and p*o — 1-2 x 10^ 
Mq/pc^. Note that p^o = 69 Mo/arcsec'^. This power- 
law distribution may persist down to the tidal disruption 
radius, Rt, which is, according to Frank & Rees (1976), 



Rt ~ 2.0 X 10'-^ cm m 



-1/3 



(15) 
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where m* is the mass of the star in Solar masses. 

For convenience we choose p = 3/2 as the resulting 
integrals are the simplest and the resulting accuracy is 
adequate for this paper. The mass of the stars within a 
given radius R is 

M{R) ~ 600 {RbRf^'^ , (16) 

where both Rb and R are in units of arcseconds. For ex- 
ample, at i? = 0.1, 1 and 10 arcsecond, we have ^ 600, 
2 X 10'' and 6 x 10^ M©, respectively. 

Number of star-disk crossings per unit time and unit 
area is n^,{R)v±{R) ~ (1/V3)7i*(i?)wx(i?) where n^{R) — 
p*{R)/m9f and v±{R) is the average of the star's velocity 
perpendicular to the disk, which we accepted to be l/\/3 
of the local Keplerian value. Note that stars with both 
positive and negative v± should be counted since both 
produce a flare. Taking the integral over the disk surface 
area from the inner radius to the outer one, the number 
of star-disk crossings per unit time is 

N{R) = (47r/V3)VGMm^n,oi?6/'ln§^ = 0.51n§^m;i 

(17) 

in units of day ^. If Rd/Ri — 500, then this estimate yields 
^ 3 star-disk crossings per day, whereas if Rd/Ri = 10 
then we have ^ 1 crossing per day. The estimate (|17ll is 
very close to the observed rate of the X-ray flares of about 
1 flare per day (Baganoff et al. 2003b). 

3.2. Ef Facts of the stellar cluster on the disk 

Under the assumption that the disk is truly inactive, i.e. 
its viscosity is zero, the only mechanism for a radial mass 
transfer in the disk is the angular momentum loss due 
to the star-disk collisions. Ostriker (1983) calculated the 
rate of the angular momentum loss by the disk assum- 
ing an isotropic star cluster. We will use his results for 
the Bahcall-Wolf (1976) distribution of stars since the lat- 
ter yields n*(i?) c>c R~'^/^, i.e. roughly same as observed 
(Genzel et al. 2003). For relatively small radii that we 
are interested in, i.e., R ^ 10*i?g, the angular momentum 
transfer from the star to the disk is dominated by the 
hydrodynamical interaction (rather than by tidal effects) 
and so we set 7y = (see equation 10 in Ostriker 1983). 
Ostriker (1983) defined the "drag time" td as the time for 
the disk at radius R to loose all of its angular momentum 
due to the star-disk collisions. For the Bahcall-Wolf (1976) 
distribution, 

0.09 

- n^AR)Rivk{R) ' ^ ' 

Note that td oc R^, and therefore it is always the inner 
disk that can be destroyed most rapidly (accreted onto 
the black hole) due to the star-disk collisions. 

As in H2.5I suppose that the age of the inactive disk is 
ta — lO^tae years. Then the disk regions with td{R) ^ ta 
will be emptied out by the star cluster and thus we can 



define the inner disk radius Ri through td{Ra) — ta'- 

- ^V o.092l%/i^. - ^ ' 

(19) 

where is the star-disk collision rate in units of 10"^ 
year~^, and where we set Ri = IQ^Rg in the slowly chang- 
ing logarithmic factor. Note that the factor in the denom- 
inator of equation H19|l is ~ 4 so that a simple interpre- 
tation of the result is that the disk area within radius Ri 
approximately equals to the area of the disk covered by 
the star-disk collisions in ta years. It is also worth observ- 
ing that the inner disk radius depends relatively weakly 
on the star-disk collision rate N and ta. 

Finally, the inner disk may be emptied by means other 
than direct stellar impacts. For example, the S2 star is as 
luminous as 10^ Solar luminosities (Ghez et al. 2003b), 
and it came to within 2 x lO'^i?^ off the black hole (Schodel 
et al. 2002). Clearly its separation from the disk surface 
would be very small at the time of the closest approach 
and the disk is heated to temperatures as much as 10^ 
K. The disk may become ionized enough to actually start 
accreting onto the black hole (at these small radii only). 
Therefore we believe that our conclusion that the inac- 
tive disk is likely to miss its inner part is quite robust. In 
addition disks in other LLAGN do have inner holes with 
R, ~ 100 to IG^Rg (Quataert et al. 1999; Ho 2003). 

3.3. Summary: the role of the Sgr A* cluster 

We have shown f tl3.1|l that the rate of star-disk crossings 
is few per day for the stellar cusp distribution given by 
Genzel et al. (2003) This calculated event rate is rather 
close to what is observed (1 flare per day). In addition, 
we found that the star impacts on the inner disk are so 
frequent for R ^ 10^i?g that the inner disk may be simply 
emptied out, i.e. accreted onto the black hole. The inner 
disk radius was estimated to be Ri ~ 10^i?g. 

4. Schematics of a star passage through the disk 

4.1. Background notes 

The subject of star-disk interactions in AGN has been ex- 
tensively studied in the literature. There were two main 
themes in these studies. One was the influence of the 
stars on the disk evolution through the additional drag 
imposed by the stars on the disk (e.g. Ostriker 1983; 
Norman & Silk 1983; Hagio 1987; see also S3, heat- 
ing (e.g. Perry & WiUiams 1993) and mass deposition due 
to stripping and ablation of star's envelopes (Armitage, 
Zurek & Davies 1996). Not less important were the stud- 
ies of the influence of the disk on the distribution of 
the stellar orbits in the star cluster (e.g. Norman & Silk 
1983; Vilkoviskij 1983; Syer et al. 1991; Ranch 1995; 
Karas & Subr 2001; Vilkoviskij & Czerny 2002). Zurek, 
Siemiginowska, & Colgate (1994) proposed that the wake 
region behind the star (when the star leaves the disk) may 
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be the putative "broad line region clouds" of AGN. Ivanov, 
Igumenshchev, & Novikov (1998) studied the passage of a 
black hole through the accretion disk. At the same time 
relatively little attention has been given to the emission 
from a single star-disk passage. The point is that this emis- 
sion is usually weak on the background of the AGN and 
is unlikely to be seen, although the effects of many star- 
disk collisions may be observable (Zentsova 1983; Perry & 
WiUiams 1993). 

Furthermore, in the midplane of the dense AGN-type 
disks, the cooling time behind the shock is very short, 
and the shocked gas reaches temperatures of "only" ^ 
10'^ — 10^ K. Thus no significant X-ray emission was ex- 
pected. However, in the case of Sgr A* the disk density 
is probably much lower than the typical values considered 
in the literature on star-disk collisions so far. In addition, 
even if the disk midplane density is high enough so that 
the disk is optically thick, density becomes quite small in 
the disk photosphere. If the density is low, the radiative 
cooling time is long, and the shocked gas is heated to tem- 
peratures as high as 10^ — 10^ K and emits X-rays. This 
emission is too weak in the context of normal AGN. For 
Sgr A* , however, the combination of the extreme dim- 
ness of its quiescent X-ray emission and its proximity to 
us makes the flares observable. 

4.2. Phases of the star passage through the disk 

The star's internal density is very much larger than that 
of the disk and so effects of the disk on the star dur- 
ing one passage are minimal (see Syer et al. 1991). An 
alternative way to see this is to compare the mass of 
the cold gas swept up by the star as it travels through 
the disk, to the stellar mass: AM/M* = J^nRl/M^ = 
1.4 X lO-^miTj/Vf V^m-i < 1. We thus consider 
the star to be a rigid solid body in this paper. We can 
also neglect accretion of gas onto the stellar surface be- 
cause the corresponding Bondi radius is much smaller 
than {Rb = GM^mp/2kTchi,^{R) = 4.5 x 10^m*r4 < 
7 X 10^°r«). Tidal interactions of the star and the gas in 
the wake of the star are hence neglected. 

The star drives a shock into the disk. The Mach num- 
ber is very large, v^/cdisk ~ v^Tkl^T^ ^ lO^i and thus 
the star drills a very narrow (i.e. with diameter ~ several 
i?*) hole in the disk. If the disk is optically thick, then the 
X-rays will be observable only when the star, or the shock 
front, is in the disk photosphere. Figure Q shows the ex- 
pected sequence of events, with time increasing from the 
bottom to the top. The star is shown as a black filled circle 
with an arrow indicating its velocity. The optically thick 
denser part of the disk is shown with a darker color, and 
the disk photosphere with a lighter one. For convenience 
we call the side of the disk impacted by the star first the 
"front side", and the other side as the "back" one. With 
a yet lighter color we show the region of the disk heated 
by compression and by radiation to temperatures of order 
~ 10^, perhaps up to ~ 10^ K. Finally, the hottest region 



Afterglow 




Thermal Flare 




Front side X-ray Flare 





T~ lO^K 


1 





Fig. 1. Cartoon of a low mass star passing through the 
disk. Star is shown as a filled black circle with an arrow. 
The disk is schematically divided into an optically thick 
disk (darker) and an optically thin - the photosphere - 
parts (lighter color). Time runs from bottom to top. The 
X-rays are observable only when the star is in the pho- 
tosphere of the disk (first and third sketches from the 
bottom). When the star is deep inside the disk (second 
sketch) , the emitted X-rays are absorbed and thermalized 
and emitted in the near infra-red. In the afterglow phase 
(top) the hot shocked gas cooled to low temperatures and 
only weak infra-red emission is present. 

is the one shown with a light colored circle around the 
star, with temperature of the order the characteristic one 
(eq.mi). 

Since the star's motion through the disk is highly su- 
personic, the disk density remains approximately constant 
during the whole event except for a small region within 
several stellar radii from the path of the star. This means 
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that the region of the disk pre-heated by X-rays and the 
shock will not have enough time to expand much and form 
a bulge next to the exit point of the star. The disk re- 
mains approximately plane-parallel during the star's pas- 
sage through it, then. A good analogy to the situation is 
given by a bullet passing through a low density material 
(e.g. a mattress). 

For the optically thick disk, the passage of the star 
through the disk can be divided into three stages (cf. 
Figure 1): 

(1) The thermal burst stage - the star is inside the 
optically thick part of the disk. X-rays are absorbed within 
the disk and re-emitted as thermal radiation in the infra- 
red or optical. This is the regime in which the star passages 
were extensively studied in the literature. 

(2) The X-ray flare phase - the star is in the disk pho- 
tosphere. The X-rays are directly visible to the observer in 
this phase. In the simple semi-analytical treatment of this 
paper, we will not make distinction between the front-side 
and the back-side flares, although in reality the observed 
emission will depend on the 3-D geometry of the problem 
through i.e. obscuration of a part of the shock by the star 
itself. 

(3) The "afterglow" phase ~ the star is many stellar 
radii away from the disk. The X-ray emitting gas cooled 
down to temperatures below ~ 10^ K and is not observable 
anymore (due to high interstellar absorption to Sgr A*). 
A patch of the disk around the trek of the star is still 
"hot" compared with the unperturbed disk, and is cooling 
mostly by black-body emission coming out in the infra-red 
(e.g. Syer et al. 1991). 

If the disk is optically thin this division does not apply 
since there is essentially only one stage. 

4.3. Division of the parameter space 

Radiation from the shocked gas is a very important facet 
of the given problem since it affects both the hydrody- 
namics and the expected spectra. Therefore, the parame- 
ter space for the problem naturally breaks down into two 
important cases: optically thin and optically thick (disk). 
We will consider these cases separately below. In addi- 
tion, there is further dependence of the results on whether 
the cooling time of the shocked gas is shorter, comparable 
or longer than the adiabatic expansion time (see below). 
Altogether we have six cases each of which may lead to 
different expected flare spectra. Some of the resulting flare 
characteristics, such as the duration, depends mostly on 
geometry rather than on anything else. 

5. X-ray flares (luminosity) 

5.1. Optically thin case 
5.1.1. Setup 

Consider the case of a disk optically thin to X-rays as 
seen by the observer. The photon energy window through 



which Chandra and XMM are able to observe Sgr A* is the 
i? ~ 2 — 8 keV range. The opacity in the disk photosphere 
is of course a very strong function of photon energy E, and 
at small energies, i.e. E ~ few kcV, photo- absorption is 
very significant and the opacity is much higher than at en- 
ergies £' ~ 7— 8 keV. We will consider this effect in greater 
detail in i j6.2l but for now we introduce an "effective" disk 
opacity, one which is a rough approximation to the opac- 
ity weighted with the optically thin (photon energy) spec- 
trum of the shock. Since we expect bremsstrahlung-like 
hard spectra (see below), and since photo- ionization of 
the gas may be important (see t i6.2|) , we estimate that the 
effective disk opacity is 6 ^ few times ctt, the Thomson 
one. 

Within this approach, the optical depth of the disk in 
the relevant photon energy range along the line of sight to 
the observer is 

r ~ —arHrid , (20) 
cos* 

which yields the maximum midplane density for which the 
disk is still optically thin: 

< not ~ 'J-86 X 10 ^ I2 ' . (21) 

Now let i?* — r, Rq be the star's radius, M* = m, M© 
is its mass and is the star's velocity. We define the 
relative velocity of the star and the gas in the disk, 
v^^Y — v'j^+v1 — 2vifVK cos Or, whcrc 9r is the angle between 
the two velocities. In general one needs to carry out cal- 
culations for arbitrary and 9^ and then integrate over 
the 3D stellar velocity distribution. Here we use the sim- 
plest approach and take the star velocity be perpendicular 
to the disk surface, i.e. Or = 7r/2. Since we also assume 

= vk = 2.1 X 10^r4 cm/sec, we have Wioi = V^vk- 
We shall emphasize that in a more detailed calculation 
the role of the exact values for Or and will be instru- 
mental in determining the fiare characteristics. Moreover, 
the exact 3D velocity of the star will also be important. 
For example, if the star is moving away from the observer, 
then the leading edge of the bow shock is obscured from 
the observer by the star. This effect is not taken into ac- 
count in our calculations. Thus we should keep in mind 
that the calculations below are meant to be representa- 
tive and in reality a larger variation in the spectral and 
temporal characteristics of flares is expected. 

5.1.2. Flare duration and rise time 

If the disk height scale is much smaller than the radius 
of the star, then duration of the burst is given by tdur — 
2_R»/u,: 

tdur - 2— = 670 sec r^r]''^ , (22) 
In the opposite case we have i^ur ~ -ff/w* 

idur - 2 — = 5 X 10^ Tl'^rl sec . (23) 
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Note that H — 2i?* at r4 = 0.1 Ir*' , so that for radii 
smaher than this one equation [53 should be used while 
for larger radii equation [53 is the one to use. Since the 
maximum luminosity reached in the bursts can be signifi- 
cantly higher than the quiescent Sgr A* emission, the rise 
time scale should be at least several times shorter than the 
burst duration. However it cannot be much shorter than 
a fraction of i.e. a hundred seconds for r^r^^ ~ 1. 



5.1.3. Characteristic temperature and time scales 

The maximum temperature to which the gas can be 
heated in the shock is the "characteristic" temperature 
defined as 



cha 



3 2k 



1.8 X IOV4 1 K, 



(24) 



where we assumed that the mean particle mass in the 
shocked gas is = mp/2. 

In a strong shock, the gas is compressed by the factor 
of 4, so that the electron gas density in the shocked gas 
is Ugh — 4nd. For the temperatures of interest the cooling 
function is dominated by the bremsstrahlung radiation. 
Accordingly, A(T) = 2.4 x IQ~'^'^T^/'^ erg cm^ sec"^ for a 
gas of cosmic abundance (we will include atomic emission 
processes in the calculation of spectra later). The cooling 
time of the shocked gas is 

= A^^T— ^ ^ ^^'^ ^4 '^'«n . (25) 
A(T)nsh 

where we put T = T'char(^)- The radiative cooling time tc 
should be compared with the adiabatic losses time scale, 
tad, which is the time it takes for the shocked gas to cool 
off via adiabatic expansion in the disk. 

iad — — — 330 sec r^r^^ . (26) 

If the gas density is very high, so that tc <C tad, then 
the shocked gas will not be heated to the maximum tem- 
perature (eq. I24|) because of the radiative cooling. We 
estimate that in this case the actual temperature of the 
shocked gas near the star is 

(27) 



cha 



^ad ^" tc 



5.1.4. X-ray luminosity 



The star is essentially a piston pushing the gas ahead (and 
aside) of it, and hence the rate at which the star makes 
work on the gas in the disk is approximately 

-3/2 



7 X 10^"* erg sec ^ n-urlr^ 



(28) 

where nn — rid/lO^^ cm This is a mechanical power, 
and only a fraction of this energy may be emitted in X- 
rays. Defining ex to be the efficiency with which mechan- 
ical power is radiated away in X-rays with E ^ 1 keV, we 
write 

Lxe = ^xLyj , (29) 



where the subscript "xe" stands for "X-ray emitted" . The 
efficiency depends on the ratio of tc and <ad- Below we 
use rather rough estimates of ; better values are unlikely 
to be obtained in an analytical way. 

If tc <SC tad, then shocked gas temperature (equation 
I27|l is smaller than Tchar due to strong radiative cooling. 
If Tsh is lower than ~ 10^ K, then most of X-ray emis- 
sion will be emitted at soft X-ray energies. Such emission 
would be absorbed in the ~ lO^'^ hydrogen nuclei per cm^ 
column depth of the cold material in the line of sight to 
Sgr A*. In other words flares with Tgh ^ 10^ K arc not ob- 
servable from Sgr A* location. However we find that the 
disk midplane density that explains the observed luminos- 
ity of flares in Sgr A* best is low enough that this situation 
(Tsh ^ 10^ K) is not reached. It is then appropriate to set 

— 1 in the radiatively efficient case. 

In the opposite limit of a very long cooling time, the 
shock becomes an inefficient radiator. The X-ray lumi- 
nosity of the shock front around the star itself is L ~ 
Lwtad/itad + tc) ~ L^t^d/tc- Thus in general we write 



tad 



^ad ~t~ tc 

In particular, in the case of tc ^ tad, we have 
Lxo 4.1 X 10^^ erg sec"^ "n?'*^^ 



(30) 



(31) 



It is important to note that equations 1301 and 1311 do 
not include contribution of the shocked gas in the wake 
of the star. When radiative cooling is initially unimpor- 
tant, the shocked gas cools first via adiabatic losses and 
work against the ambient unshocked gas. As the hot gas 
expands, the adiabatic cooling time scale tad eventually 
becomes comparable to the radiative cooling time tc, at 
which time most of the radiation would be emitted. If 
the gas temperature at that point is still larger than 10^ 
K, then the X-ray luminosity could be as large as L^,. 
Unfortunately we found no reliable analytical way to accu- 
rately calculate this effect and we leave this task to future 
work, remembering the fact that tx is somewhat underes- 
timated in the tc S> tad case. 

5.2. Optically thick disk case 

We now assume that the disk optical depth to X-rays as 
seen by the observer is greater than unity, and hence (cf. 
equation [2J 



nd>not = 3.86xl0"^r-'^V^' 



(32) 



The main difference from the optically thin case is the fact 
that X-rays emitted by the shock while the star is in the 
disk midplane are absorbed and reprocessed into thermal 
radiation. As explained in H4.ll this is the regime studied 
intensively in the previous literature. However even in this 
case there will be X-ray emission when the star reaches the 
disk photosphere. There the density is usually low enough 
to preclude radiative cooling from reducing shock temper- 
ature to below 10^ K. Hence the ratio of the fraction of 
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energy emitted in X-rays to that in the optical-UV in the 
optically thick case is 



fdtL^ 1 

J < - « 1 , 

J at Louv Td 



(33) 



where Td is the optical depth of the disk to X-rays. Note 
also that two X-ray flares are emitted per each star-disk 
encounter - one at the impact side and the other at the 
exit side of the disk (see Fig.QJ, but of course the observer 
will see only the flare emitted from the side of the disk 
facing the observer. 

The flare duration is now at most H/v^, because the 
flare from only one side is seen. In addition, the duration 
will be reduced by a further factor of a few if the photo- 
sphere of the disk is at several hight scales above the disk 
(since the disk density drops very quickly there; see eq. 



Lxo — ^xc oc exn[z) (equation I29|) . Later, when the sur- 
face Te = 1 is reached, the attenuation of X-rays by the 
disk matter becomes important and the luminosity de- 
pendence on the density becomes Lxo oc txn{z) eyij>{—Te) . 
Furthermore, for the gaussian density profile 



Te{z) 



n{z) 
rid 



(36) 



where t is the optical depth to the disk midplane (equation 
I2UII and is the disk midplane density. It then follows 
that 

/ \ / \ r "-(^) 

,[t) (X €x n[z) exp ~ 



rid 



(37) 



Thus the event is indeed a flare: the luminosity first in- 
creases and then decreases quite rapidly. Note however 
that we implicitly assumed here that cooling time tc is 
shorter than duration of X-ray flare, tdm- 



5.2.1. X-ray light curve 

The emitted X-ray luminosity can be calculated exactly 
as in the optically thin case f H5.1.4|l but we now take into 
account (i) attenuation of X-rays in the disk; and (ii) the 
fact that unshocked gas density in the disk photosphere, 
rip, is lower than that in the disk midplane, rid. Our treat- 
ment of the first complication is to simply introduce expo- 
nential attenuation along the line of sight. This approach 
underestimates the X-ray luminosity for very small cos i 
for which radiation scattered once in the photosphere be- 
comes more important than the unscattered one. We defer 
a more careful treatment of the scattered radiation to a 
future paper however, to keep our calculations here as an- 
alytical and transparent as possible. 

One should also keep in mind that we use the "effec- 
tive" cross section for the total one (photo-absorption plus 
scattering), which we parameterize as bcTT (see H5. 1.1(1 . 
Now let Te be the effective total depth of the disk from 
the star's position z to the observer. The observed X-ray 
luminosity is then 



Lxo ^ exp [-Te] , 



(34) 



where Lxo is given by equation 1291 The gas density in 
the photosphere is approximately given by the gaussian 
density profile: 



n(z) = rid exp 



H 



(35) 



The observed X-ray luminosity is now a function of 
time since z = v^,zt if we define t = as the time at which 
the star crosses the disk midplane. Consider a front side 
flare - when the star is moving into the disk (i.e., see the 
sketch on the very bottom of Fig. ^ the observer is be- 
low the disk, at z = — oo in this case). When the star 
first impacts the disk photosphere, the X-ray luminos- 
ity rapidly increases because the star moves into denser 
and denser environment. During this optically thin stage 



5.3. Maximum flare X-ray luminosity 

The maximum observed X-ray luminosity, Lmax, for an op- 
tically thick disk, will be approximately that correspond- 
ing to Te = 1. The gas density at which Te = 1 is simply 
n{z) = riot as given by equation (|^ . We now substitute 
this value for the gas density into equation[2niand multiply 
the result by exp[— 1] to get 



10^" ex(not) 



cos I -1/2 2^-3 

— J 2 r,r^ 



(38) 



where b ~ few is the ratio of the effective total cross sec- 
tion in the disk atmosphere to the Thomson cross sec- 
tion. Note that in the case of an optically thin disk the 
midplane density is lower than riot and hence the emit- 
ted X-ray luminosity is lower. Therefore eauationl38l gives 
the maximum X-ray luminosity of a flare in both optically 
thin and thick cases. 

We shall emphasize the fact that the maximum X-ray 
luminosity just derived is a strong function of the disk in- 
clination angle i. In particular in the radiatively inefficient 
shock case, oc cosi and thus L^ax oc cos^ i. This depen- 
dence shows that if X-ray flares are indeed associated with 
the putative inactive disk, then the inclination angle could 
not be too large, i.e. the disk could not be face on or else 
the X-ray luminosity of the flares would be too small. 

6. X-ray spectrum 

6.1. Optically thin spectrum 

In general the X-ray spectrum of a flare is complicated. 
The emitting region is a multi-temperature gas and only 
numerical calculations are likely to yield accurate spec- 
tral predictions. However single temperature optically thin 
spectra will provide us with certain guidance. In addition, 
these spectra should be of some relevance for the optically 
thick disk and the case of tc ~ tad • Under these conditions 
we expect Tgh ~ LLhar and a major fraction of the me- 
chanical power of the shock should be "quickly" radiated 
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in X-rays, so that the wake of the star does not make a 
large contribution to the emitted spectrum. 

The photo-ionization parameter for the gas in the 
shock is ^ ~ Lx/rishB? and is usually smaller than 
a hundred, i.e., photo-ionization is not important for 
the shocked gas, and the plasma is essentially in the 
coronal limit. The ionization equilibrium assumption is 
appropriate since the recombination time scale for the 
completely ionized Fe is l/(arec'^sh) ^ 3 sec n'^lT^^'^ . 
Further, the Thomson optical depth of the shocked gas 
is tt ~ arnshR* = 0.02r*nii <C 1, so that the emis- 
sion region is optically thin. The Compton y-parameter, 
y = A{kTsh_/meC^)TT <C 1 in most realistic cases and thus 
Comptonization is not important for the continuum. 

We computed optically thin X-ray spectra using the 
code XSTAR (e.g., Kallman & Krolik 1986), assuming so- 
lar abundance of elements and that the emitting gas tem- 
perature is equal to Tchar(-R)- The results for three differ- 
ent radii, — 0.5, 1 & 2.5, are plotted in Figure |21 with 
the dotted, thick and thin solid curves, respectively. The 
near infrared part of the spectra will be discussed later 
in 13 In Figure 12 the top panel shows the spectra in a 
broad frequency range and the bottom panel concentrates 
on the 1 — 10 keV part of the spectrum. 

The characteristic temperatures are Tchar = 
3.6,1.8 and 0.73 x 10® K for n = 0.5,1 and 2.5, 
respectively. For such high temperatures, the X-ray 
spectrum is mostly bremsstrahlung emission with the 
great exception of the Fe Ka lines at ~ 6.7 and 6.9 keV. 
The keV-part of the spectrum shows that the line is 
dominated by the He-like (i.e., the 6.7 keV) component 
for the smallest T considered, whereas the 6.9 keV 
component dominates for the largest value of T . 

Clearly, the integrated equivalent width (EW) of 
Fe Ka line decreases as Tchar increases (i? decreases). 
Another obvious trend seen in Figure |(5J) is that with in- 
crease in Tchar the thermal rollover moves to higher fre- 
quencies and the spectrum begins to look more and more 
power-law like in the 2 — 8 keV band. To quantify these 
trends and to compare the observed flare spectra with 
those predicted by our model, we defined the photon index 
r by drawing a power-law through two line- free photon en- 
ergies oi E ~ 2.2 keV and 7.25 keV, and also integrated 
over all the components of the Fe Ka line to get its total 
EW. Figure 13 displays the resulting dependence of index 
r and the EW of Fe Ka line on the temperature of the 
shocked gas. 

So far only two - the brightest - of the observed X- 
ray flares had enough statistics to conclude that spectra 
during flares seem to be hard power-laws in the 2 ;S ^ 8 
keV energy range, with the photon spectral index F ~ 
1±0.7 (Baganoff et al. 2001, Goldwurm et al. 2003a), and 
they seem to show little of an Fe Ka line. Goldwurm et 
al. (2003b) notes that a Fe Ka line-like feature at ~ 6.4 
keV in their spectra is not statistically significant; but the 
upper limit on the EW of the line is as large as 1.8 keV. 
From Figures |21 & El it is clear that to have such hard 
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Fig. 2. Optically thin X-ray spectra for three different val- 
ues of disk radius: — 0.5, 1 and 2.5 for the dotted, 
solid thick and thin curves, respectively. The gas density 
is rid = 10^^ cm^'^, inclination angle i = 70°, and the star's 
radius is i?* = 2 Rq for all the three cases, (a) - Spectra in 
a broad frequency range. The near infra-red curves will be 
explained later in fJ7| (b) - the X-ray part of the spectrum. 
The upper curve is shifted up by 3 whereas the lower one 
shifted down by the same factor for clarity. Note that, ex- 
cept for the brightest of the curves, the "hot" Fe Ka line 
is clearly present. 

spectra, the bright flares should have Tgh ^ 2 x 10® K. 
The constraint on the EW of the Fe Ka line also requires 
similarly high values of Tgh- 

Finally, let us recall that in the case tc 3> iad, we 
expect that there will be a range of temperatures that 
contribute to the overall spectrum, from T = Tchar in front 
of the star to much lower temperatures in the wake of the 
star. The main effect is to make the spectrum softer. Hence 
the minimum value of F ~ 1.3 from Figure El should also 
be the minimum value for the spectral index in the more 
general case. 

6.2. Absorption along the line of sight in the 
photosphere 

The emitted spectrum is attenuated in the atmosphere 
differentially for different photon energies E because pho- 
ton opacity is a function of E. The observed spectrum 
could be severely cut (compared to that shown in Fig. 



12 



Nayakshin et al.: X-ray Flares and inactive disk in Sgr A' 



Opacity (Thomson units) 



Radius, Rg 
,4 




Fig. 3. (a) log(EW) of Fe Ka line and the X-ray spectral 
index F versus temperature Tsh of the emitting gas for 
optically thin X-ray spectra for Solar abundance of ele- 
ments. If Tsh = 7char(-R), then F and EW may be plotted 
versus radius R, shown as the upper x-axis. Note that to 
yield Fe Ka EW and spectral index F compatible with the 
observations of the two bright flares (Baganoff et al. 2001 
& Goldwurm et al. 2003a), T^h ^ 2 x 10* KotR< 10* Eg. 

12) at low energies E ^ 2 keV because the opacity there 
is very large. Observations of at least the brightest flares 
show no enhancement in the column depth of the cold 
neutral absorber in the line of sight (Baganoff et al. 2001, 
Goldwurm et al. 2003a). For optically thick case, the only 
way to reconcile the model with the observations is if the 
photosphere around the star is strongly ionized and hence 
the absorption opacity is much smaller than that due to 
the cold material in the line of sight. 

Baganoff et al. (2001) found that the column depth of 
the cold material during both quiescent and flare emission 
is iVcoid 5 X 10^^ hydrogen nuclei per cm^. The column 
density through which the X-rays pass during the maxi- 
mum X-ray observed emission is around A'ion ~ l/^c* ~ 
5 X 10^'^. We then need the absorption opacity of the ion- 
ized photosphere to satisfy 



0.1 



(39) 



where (Xcoid is the absorption cross section of the cold gas 
which is much larger than ax, the Thomson cross section. 
For energy E in equation 1391 it is sensible to take i? ~ 2 
keV because the cold gas opacity is largest at the smallest 
energies in the keV band, and so the measurement of A^coid 
is really dominated by the lowest energy bin (see Figure 3 
in Baganoff et al. 2001). 
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Fig. 4. Photo-absorption cross section (in units of ctt) 
of photo-ionized matter as a function of photon energy. 
Values of photo-ionization parameter ^ are indicated next 
to the respective curves. Note a very fast decrease in 
the opacity with increase in the ionization parameter for 
^ ^ 1. This occurs because of photo-ionization of the ele- 
ments providing the opacity. 

The ionization state of the disk atmosphere illumi- 
nated by X-rays from the shock is controlled by photo- 
ionization rather than by coUisional ionization since the 
gas temperature is low, i.e. T ~ T^r ^ lO'^ — 10'' Kelvin. 
Using XSTAR once again, we calculated the ionization 
state of plasma kept at T = 10^ K and illuminated by 
X-rays with spectrum similar to those shown in Figure 
(0). The resulting absorption cross section as a function 
of photon energy is plotted in Figure Q for different 
values of photo-ionization parameter which is defined as 
^ = An Fx /no where F^ is the integrated X-ray flux. The 
solid curve labelled "0" is calculated for a very small ^ 
and yields the "cold" absorption cross section crcoid(-E'). 
The Figure shows that the condition specified by equa- 
tion H39|) is fulfilled when ionization parameter exceeds 
about 10. 

We estimate the average ionization parameter in the 
photosphere around the star as 



no(i//2)2 



(40) 



where Lxe is the emitted X-ray luminosity at the max- 
imum light, which is given by Lmax (see SQ, and as 
representative distance from the star H/2 is chosen. For 
example, Lmax for the case of an optically thick disk with 
adiabatic losses dominating the radiative ones is 



0.042^r-^/Vfr7^ 



(41) 



Clearly strong enough ionization of ^ ^ 10 requires the 
two bright detected X-ray flares to happen at small values 
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of r4 and/or be produced by large stars with ^ 1. 
Taking a concrete example of a 5 Solar radii star, we find 
that ^ exceeds 10 for radii r4 ^ 0.6. Note that the photo- 
ionization parameter at the time of maximum light in 
X-rays is a function of the disk inclination angle, so once 
again we see the dependence of the results on the disk 
orientation to us (cf. ^5.2.1l and l5.1.1|) . 

We shall emphasize that the constraint ^ ^ 10 is im- 
portant only if the disk is somewhat optically thick. The 
disk may be optically thin to X-rays which we find quite 
likely for the innermost disk region where the bright flares 
should be taking place. 

7. Flares in the near infrared 

Here, as everywhere else in this paper, we concentrate on 
the effects of relatively low mass stars on the disk. For 
high mass, high luminosity stars, the main effect through 
which the disk is affected is the heating of the disk by 
the enormous optical and UV radiation output of these 
stars (Cuadra et al. 2003). In addition winds from such 
stars may increase the effective radius of the star: due to 
the winds a larger area of the disk may be shock heated. 
These stars should be expected to produce NIR and X-ray 
flares much stronger than those we consider in this paper. 
Of course such powerful flares are much less frequent than 
those from low mass stars, and we plan to consider in our 
future work. 



7.1. Optically thick case 

First consider the simplest case of a disk optically thick in 
the near infrared. X-rays emitted by the shock while inside 
the optically thick disk are easily absorbed by the cold 
un-shocked gas and reprocessed into thermal radiation. 
Therefore we assume that all the work done by the star 
on the cold gas is emitted as a blackbody radiation in this 
case and that the luminosity is equal to ~ (equation 
I28() . The radiation runs ahead of the shock wave, heating 
and ionizing the gas in the disk. The region affected by 
the radiation is roughly the cube with size 2H. The 
temperature of the region can be found by assuming that 
the cooling time within the disk, ted, is shorter than H/v^. 
This imphes that ~ ctbT^, x 2{2Hf and 



T,r = 1.78 X lO-' K 



-,1/4 



-9/8 



(42) 



(index ir stands for "infra-red"). If t^d > t-ptc, then the 
thermal luminosity is lower than L^. 

As remarked in ^ the pre-heated gas cannot expand 
much during the star's passage through the disk. The ob- 
served flare luminosity in near infra-red is then reduced by 
the disk projection effect: iobs = cosi. The near infra- 
red K-band seems to yield the tightest constraints on the 
theoretical models (see, e.g., Hornstein et al. 2002). The 
corresponding reference wavelength is 2.2/im and the fre- 
quency is — 1.36 X 10^** Hz. The predicted observed 



luminosity is then vLi, = 47r x AH^ cosz vBi,(Tir): 

36 . 



vL^ = 7.8 X lO'^'cosi T2 



- 1 ' 



(43) 



where v = v/vq and x = hv/kTir, and where we as- 
sumed that the distance to the Galactic Center is D = 8.0 
kpc. Note that the factor of 4tt above arises due to the 
usual assumption of the observed at Earth that the emit- 
ter (the disk) emits isotropically in the full An steradi- 
ans. Also note that the dependence on the disk density is 
"hidden" in the temperature T^j, of the emission and that 
x = i? {hvo/kT„) = v (6.5 X WK/T„). 

The luminosity given by eauation l43l is quite large. The 
current best limit on the quiescent emission of Sgr A* at 
2.2 /im is about 2 mJy (deredenned; Hornstein et al. 2002). 
In terms of vL^,^ the quiescent limit of 2 mJy corresponds 
to about 2x 10"^"' erg/sec. Hornstein et al. (2002) note that 
the limits on random 3-hour flares is about a factor of 10 
higher, i.e. ~ 2 x 10'^^ erg/sec. Hence an optically thick 
(in NIR) disk in Sgr A* would produce NIR flares strong 
enough to be detected by now. 

7.2. Optically thin disk 

Previously we pointed out that no optically thick disk can 
exist in Sgr A* already based on its quiescent spectrum - 
or the disk would reprocess too large a fraction of the 
visible and UV stellar flux into the infrared band and 
would violate the observed upper limits (see ti2.2|) . These 
constraints limit the infrared optical depth of the disk to 
Tir ^ 0.01. If the opacity of the disk grains were grey, the 
grain temperature would equal to the effective one (that 
is Tir). In this approximation the luminosity of the flare 
in the NIR is simply reduced by the factor Tii-: 



7.8 X 10 



34 



0.01 



cos i T2 r| 



1 



(44) 



Figure 121 panel (a), shows examples of such spectra for a 
star with radius R — 2Rq. Note that the NIR flares are 
weak enough to be missed. Also note that the value of 
the dust optical depth, assumed to be equal to 0.01 for all 
frequencies, is likely to be smaller than that of the disk in 
quiescence. The stellar passage through the disk may in 
fact destroy the dust (see below) in which case the NIR 
luminosity will be reduced even further. 

It is tempting at this point to make a connection be- 
tween S, the disk surface density given by equation El 
and the optical depth in the NIR using the standard in- 
terstellar dust opacity. For example, at 2.2/Lim, the dust 
extinction is about 5 x 10~^^ cm^ per Hydrogen atom (see 
Fig. 2 in Voshchinnikov 2003). Thus with na ~ 10^^ cm~^ 
the disk would be optically thick in the 2.2/zm for radii 
greater than about 10'^ Rg. However the dust opacity in 
the disk is likely to be reduced compared with the inter- 
stellar values but it is hard to calculate this process model- 
independently. One uncertain quantity is the minimum 
size of the dust grains which is influenced by dust growth. 
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For example, following Spitzer (1978), we estimate that 
the grain radius a grows at the rate da/dt ~ 10~'^^Qnj']^ 
cm/year, where ^ ^ is the average fraction of atoms 
that stick to the grain in the course of collisions with the 
grain. This rate is very large so there may be no small 
grains in the disk which would argue for an optically thin 
disk. On the other hand dust grains will be destroyed or 
broken up into smaller grains by dust-dust collisions. 

Another uncertainty is the ratio of the mass of the 
matter in the grains to that in the gas. There could be 
no dust in the disk initially when the active phase turned 
off because the disk was too hot; the dust could then be 
created in the disk and could also be brought in by the 
stellar winds; it could also be destroyed because of stellar 
UV radiation, star passages. For an example we briefly 
consider the last effect mentioned above. The tempera- 
ture of the gas in the cylinder with radius H around the 
star's trek in the disk is given by equation 1421 This tem- 
perature is high enough to destroy all or most of the dust 
species. Therefore we may define time "dust destruction" 
time scale, idd, as time needed to cover all of the disk 
area within radius R by star-disk collisions (each collision 
"covers" area ~ nH^ of the disk). Using the rate of the 
star-disk passages calculated in il3.1l we obtain 

idd = -■ ~ 2 X 10" years T.^Vr^ . (45) 

This time is very short compared with the viscous time 
scales for the cold disk (~ 10'' — 10* years). 

8. Typical "weak" and typical "strong" flares 

Given the theoretical picture of a star-disk collision devel- 
oped in this paper, we now consider two important cases. 
The first one is the passage of a "low" -mass star with 
i?* = i?0 through the disk for a range of midplane disk 
densities rid, and for diff'erent radii R. The second case is 
that of a "large" star, with R = Si?©, passing through 
the same disk. 

The low-mass star case is analyzed in Figure The 
midplane disk density rid is chosen to be rid = 3 x 10^°, 
10^^ and 3 x 10^"'^ hydrogen nuclei per cm'^. The X-ray 
luminosities are shown with thick curves while the near 
infrared 2.2/i luminosities are plotted with thin curves in 
panel (a). Larger values of these luminosities correspond 
to larger values of the midplane density. Panel (b) of the 
Figure shows three time scales - the duration (eq. 123(1 . 
cooling (eq. I25|l and adiabatic expansion time (eq. 126(1 . 
Out of these only the cooling time depends on the disk 
density, which was chosen to be rid = 10^^ cm~^ for the 
respective curve in panel (b). We do not show disk radii 
smaller than lO'^i?^ because it is rather unlikely that the 
disk will extend there (see tl3.2(l . 

The X-ray luminosity curves have a common asymp- 
totic at large radii, which is simply the optically thick case 
given by equation l(38() . It is very important to note that 
in this optically thick limit the X-ray luminosity is an ex- 
tremely strong function of radius: Lxo oc R~^ . Given the 
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Fig. 5. Flare properties for i?* — Rq, disk inclination an- 
gle of i = 75° and b = 2 (see eq. 138(1 . (a) - Maximum 
X-ray (thick) and NIR (thin) luminosity during flares as 
a function of radius for three different values of disk mid- 
plane densities, rid = 3 x 10^°, 10^^ and 3 x 10^^ hydrogen 
nuclei per cm'^. For NIR luminosities, we assume that the 
dust optical depth is 0.01, as in equation El (h) - Burst 
duration (thick solid), adiabatic expansion time (dotted) 
and the cooling time (thin solid). Only the latter depends 
on the disk density, which was chosen to be rid — 10^^ for 
the thin solid curve. 



fact that quiescent X-ray luminosity of Sgr A* is about 
Lq ^ 10^'^ erg/sec, it is clear that flares from solar-sized 
stars will not be observable for R >, 10^ Rg. A comparison 
with the case of a large star. Figure El shows that flares 
from these are observable up to radu of ^ few x W^Rg, 
i.e. just a factor of few higher. Moreover, this conclusion 
is independent of the disk midplane density: the X-ray lu- 
minosity of stars passing through the disk at a given radii 
cannot exceed imax- Thus we conclude that as far as X- 
ray flares are concerned, only the "inner" disk region with 
R < few X 10'' i?g matters. 

At smaller disk radii, the disk becomes optically thin 
to X-rays along the line of sight to the observer and 
at this point the X-ray luminosity curves change their 
slope. The X-ray luminosity dependence on the radius is 
much weaker in the optically thin regime (sec eq. I31() . 
Clearly "large" stars produce higher luminosities than do 
the "small" ones, in both optically thin and thick regimes. 
This is easily understood from the fact that the total work 
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Fig. 6. Same as Figure but for a high mass star with 
R = 5 i?0. Note the change in the luminosity scale; X-ray 
flares are roughly 100 times brighter than they were in 
Figure El The NIR emission increased by a much smaller 
factor mostly because it is in the Rayleigh- Jeans regime. 

done by the star on the disk is proportional to the area of 
the star feauation I28|l . plus radiative cooling is actually 
more important for larger stars (equation 13 1() . There is of 
course fewer high mass stars than there is low mass ones 
in Sgr A* star cluster (Genzel et al. 2003), and this means 
that the prediction of the model is that high luminosity 
flares should be a minority. 

Spectra emitted by flares from passages of low mass 
stars are likely to be relatively soft, i.e. with F 2. We 
reach this conclusion based on the fact that radiative cool- 
ing time for such flares is always much longer than the adi- 
abatic expansion time (see thin solid versus dotted curves 
in Fig. ISJ, and hence the gas will expand and cool sig- 
nificantly before a large fraction of the shock mechanical 
power is radiated away. As a result one should also expect 
a strong "hot" Fe Ka line to be emitted by such flares. 

On the other hand, for high mass stars the adiabatic 
expansion time is longer and hence the "average" tem- 
perature of X-ray emitting gas should be larger than that 
from the smaller stars. This is especially true for star-disk 
collisions at radii ^ few x lO'^i?^ because then the char- 
acteristic temperature is very high, i.e. Tchar ~ 10^ K, 
and even if the emitting gas temperature is significantly 
smaller than this it may still be high enough that the 
2-8 keV part of the spectrum is a pure power-law with 



r ~ 1.5 — 2. Correspondingly Fe Ka line emission should 
be much weaker from such flares. Since these are also the 
brightest flares we conclude that, statistically speaking, 
bright flares should have harder spectra than the less lu- 
minous ones, and they should show less of the Fe Ka line 
emission. 

Duration of typical flares is in the range of one to tens 
of kilo seconds. We should stress that this result is quite 
model independent. Flares with longer durations are not 
detectable because they are too weak; flares shorter that 
about a thousand seconds may not exist because we found 
that the disk cannot extend much into the innermost re- 
gion (see t|3.1|l . Even if it did, the cooling time there be- 
comes very long and the X-ray emission would then con- 
tinue for times longer than tdm- ^ -ff/w*. 

The near infrared (2.2/i here) emission is quite weak in 
all the models, at or below the '--^ 2 x 10'^^ erg/sec upper 
limit set by Hornstein et al. (2002). This is partially be- 
cause the black-body emission is a rather narrow feature 
and in most cases the 2.2/i is either on the Wien exponen- 
tial tail or on the Rayleigh- Jeans part of the curve. We 
should also note that the NIR luminosities are calculated 
under the optically thick assumption which may be vio- 
lated for the smallest values of the disk surface density 
S. However most likely optically thin spectra would peak 
in the optical band again producing little NIR flux to be 
detected (see Next "complication" in the way of de- 
tecting NIR flares from the GC in the context of our model 
is the fact that they happen not exactly at the Sgr A* po- 
sition. The angular separation correspond to i? = 10^i?g 
is about 0.1". With current telescopes a star's position 
can be measured with a better accuracy, e.g. to within 
about 0.01"(Sch6del et al. 2002), and hence if a flare in- 
deed happened ~ 0.1 arcsec away from Sgr A* position 
it would not be even counted as the one coming from the 
black hole and probably be discounted as noise. 

The cyclotron frequency in the hot blob is about 1 
GHz, assuming equipartition magnetic field. However, its 
intensity is many orders of magnitude larger than the 
Rayleigh-Jeans intensity at this frequency, so this emission 
is completely self-absorbed and therefore impossible to ob- 
serve. It is unlikely that non-thermal mechanisms could 
significantly alter this conclusion. Therefore our model 
predicts no radio-flares correlated with X-ray flares (but 
star-disk interactions may have indirect implications for 
the jet and hence some radio variability may in fact be 
due to stars; see also Nayakshin & Sunyaev 2003). 

9. Comparison with observations of X-ray flares 

Observational facts for X-ray flares from Sgr A* may be 
summarized as following (Baganoff et al. 2001, 2003b; 
Goldwurm et al. 2003a). The maximum X-ray luminos- 
ity detected so far is ^ 10^^ erg/sec. Such strong flares 
are relatively rare; i.e. there is probably ~ 5 — 10 times 
more of the weaker ^ 10'^^ erg/sec flares. Counting all 
detectable flares, there is about one flare per day on av- 
erage. There seems to be no NIR or radio counterparts to 
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the flares (BaganofF et al. 2003b; Hornstein et al. 2002). In 
particular the radio source has never varied by more than 
a factor of 2 during many years of observations (Zhao, 
Bower, & Goss 2001). The NIR limits on random ^ 3 
hour flares are about ^ 2 x 10'^^ erg/sec, much too tight 
for many existing models of Sgr A* (Hornstein et al. 2002). 
The spectra appear to be harder than those during qui- 
escence and show little if any Fe Ka line emission (both 
Baganoff et al. 2001 and Goldwurm et al. 2003a, b X-ray 
spectra contain weak emission-line like features at ener- 
gies not too far from the expected Fe Ka line but detailed 
modeling is difficult due to poor statistic. Goldwurm et al. 
2003b put an upper limit of 1.8 keV on a 6.4 keV Fe Ka 
line feature in their spectrum of an X-ray flare). In addi- 
tion there are indications that some flares may be much 
softer and are complex around the Fe Ka line energy (P. 
Predehl, private communication). To this list of Sgr A* 
flare properties we should also add the fact that if Sgr A* 
is indeed a Low Luminosity AGN, then it is the only one 
that shows such large magnitude flares. 

The following flare characteristics appear to have a 
reasonable explanation in our model (we place items in 
order from very model independent to less so): 

— (i) Frequency of flares, i.e. number of flares per day, 
which in our model is around few per day (see ^'6.1\ 
and in particular eauation ll7|l . 

— (ii) Duration of X-ray flares is in the range of ~ one 
to tens of kilo seconds. 

— (iii) X-ray luminosities of observed flares in the peak 
of their light-curves are consistent with the maximum 
allowed luminosities in our model (equation 1381 with 
ex given by equation I3U|) . Namely, the observed lumi- 
nosities are best explained by an optically thin disk 
with midplane density of order 10^^ hydrogen nuclei 
per cm^. 

— (iv) X-ray flares from similar star-disk encounters in 
the nuclei of other LLAGN. On the one hand, the max- 
imum X-ray luminosity of flares at a given dimension- 
less distance r4 from the black hole is that given by 
equation 1381 which behaves as Lx oc n^^ oc cx 
1/M|jj. That is more massive black holes produce 
weaker X-ray flares. Further, the background quiescent 
X-ray luminosity of confirmed LLAGN is always larger 
than that of Sgr A*, usually by some 3-6 orders of 
magnitude. On the other hand, some nearby Galaxies, 
e.g. M31 is not detected in X-rays, whereas M32 has 
just been detected with Lx ~ 10^^ erg/sec (Ho et al. 
2003). Hence flares in the nearby galaxies could be de- 
tectable if similar inactive disks exist there. Clearly 
the prospects of such discoveries in the future, in par- 
ticular with XEUS and Constellation-X missions, are 
much greater than they are with Chandra now. 

— (v) X-ray spectra are expected to have a distribution 
of properties. The brightest flares however should be 
produced by large stars striking the disk at relatively 
small radii, and these should have spectra as hard as 
r ~ 1.5 with little Fe Ka line emission. Less luminous 



flares may have flatter spectra, F ^ 2, and stronger Fe 
Ka line emission. 

— (vi) Near infrared (NIR) luminosity of flares produced 
by the shock around the star is around or smaller than 
the current detection limit at 2.2fi (see ^7^. 

— (vii) No flare radio emission is expected in our model 
during the star-disk impact. Such emission is com- 
pletely self-absorbed in our model. 



10. Shortcomings and limitations 

This paper presented a rather simplified treatment of the 
problem. While aware of many of the potential compli- 
cations or modifications to our results, we kept the treat- 
ment at the present level of complexity to have a relatively 
simple but broad view of the different sides of the model. 
Nevertheless it is important to point out some of these 
effects that we plan to study in the future. 

First of all, as mentioned in US. 1.11 the exact geometry 
of the star-disk passage is crucial for quantitatively accu- 
rate results. Namely, the angle that the star makes with 
the disk rotation velocity, 9^, controls the total mechanical 
energy deposition rate. Here we simply assumed 9r = tt/2. 
Further, the angle that the star makes with the normal to 
the disk determines the velocity component perpendicular 
to the disk, which is important for the duration of the X- 
ray burst (e.g. see equation 1231) . In addition, the relative 
position of the observer to the disk and the star's direction 
of motion is also relevant. The emission of the shocked gas 
viewed from the back side of the star is different than that 
from the front side of the star. 

Next complication is due to the fact that stars are 
of course not on circular Keplerian orbits. In fact many 
of the observed stars have highly elliptical orbits, many 
with e ~ 0.9 (e.g. Genzel et al. 2003; Ghez et al. 2003b). 
Therefore one needs to also allow for variation in the ab- 
solute value of u» instead of simply using = vk as we 
have done in this paper. Similarly, we did not consider the 
case of a very bright star such as S2 (e.g. Schodel et al. 
2002), or stars that produce substantial winds. Another 
complexity is that the abundance of the putative disk in 
Sgr A* need not be Solar. In fact if the disk is currently 
fed by the hot Helium rich stars, then the Hydrogen in the 
disk may be under- abundant. Additionally, Sgr A* star 
cluster may house a number of compact objects - white 
dwarfs, neutron stars and stellar mass black holes. These 
would produce flares quite different from those that we 
have studied here. 

The aforementioned effects clearly demonstrate that 
our results may be hoped to be representative of only the 
most frequent cases. Due to these effects the spread in the 
properties of the expected flares should be expected to be 
quite significant. These issues should be addressed in the 
future with a more detailed likely numerical model. 
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11. Discussion 

In this paper we proposed that the observed X-ray flares 
in Sgr A* are emitted by bow shocks around stars when 
the latter pass through an inactive disk. Such a disk has 
recently been invoked by Nayakshin (2003) to explain the 
quiescent spectrum of Sgr A* . The disk may be a remnant 
of the past accretion and star formation activity in our 
Galactic Center ( H2.ll see also Levin & Beloborodov 2003). 
Due to quiescent Sgr A* spectral constraints, the disk tem- 
perature is restricted to be of order 10^ K, whereas the size 
of the disk is at least few x 10"* gravitational radii. The disk 
should be optically thin in mid infrared frequencies or else 
it would violate the existing limits on the quiescent Sgr A* 
luminosity (see H2.2I Cuadra et al. 2003). From modeling 
X-ray flare luminosities we concluded that the disk mid- 
plane density is of order 10^^ Hydrogen nuclei per cm"^. 
Such a disk is quite "light" - it is not a subject to gravita- 
tional instability and it will not gravitationally influence 
the dynamics of Sgr A* stellar cluster significantly (see 

Using the latest results on star content in Sgr A* stellar 
cluster (Genzel et al. 2003), we calculated the rate of the 
star-disk crossings to be few per day (Q. We then found 
that the rate of star-disk crossings in the inner disk region 
with size Ri ~ 10'^i?g is so large that the disk is very likely 
to be completely emptied out (e.g. accreted) by the star 
passages within astrophysically reasonable time of ^ 10^ 
years. 

In ^ we have shown schematics of a single star-disk 
passage, and pointed out the relation of our work to the 
previous extensive literature on the star-disk passages. We 
then studied properties of the expected X-ray flares in ^ 
first in the optically thin regime f H5.1|) . and then in the 
regime in which the disk is optically thick to X-rays. In 
the optically thin disk case most of the emission during 
the star-disk passage occurs in the X-ray band. In the 
optically thick case, on the other hand, no X-rays are ob- 
servable from the shock front while it is inside the disk. 
The emission is then reprocessed in the infrared. X-ray 
flares are nevertheless are emitted when the shock arrives 
at the disk photosphere. The maximum X-ray luminosity 
is emitted when the shock front is at t ^ 1 where t is 
effective optical depth to X-rays along the line of sight. 
This then sets a very important limit on the maximum 
X-ray luminosity that could be emitted in the star-disk 
passages (equation I^HJ . 

The predicted spectra range from very hard 
(bremsstrahlung dominated) spectra with spectral index 
r ~ 1.5 to softer ones for weaker X-ray flares with F ^ 2 
( H6.1|) . The hardest spectra are also predicted to show lit- 
tle Fe Ka line emission whereas softer flares should gen- 
erally show the Helium and Hydrogen like components 
of such lines. Finally, these lines could be significantly 
red- or blue-shifted depending on whether the flare occurs 
on the approaching or the the receding side of the disk: 

— 1/2 

ASiinc ~ (vk/c) sinz 6.7 keV = 50 sini eVr^ . This in 



principle can be observable with XMM and in future with 
Astro-E(2). 

Comparing this maximum luminosity to the observed 
X-ray flare luminosities we concluded that the disk must 
be somewhat optically thin to X-rays at least in its in- 
nermost region. This then explains why there is no ex- 
cess absorption in the cold disk along the line of site in 
our model, in accordance with the observations (Baganoff 
et al. 2001, 2003b; Goldwurm et al. 2003a). In addition, 
we have shown that the disk unshocked material is also 
photo-ionized by X-rays from the shock f H6.2(l . For bright- 
est {Lx ^ 10'^^ erg/sec) flares the photo-ionization alone 
could reduce the absorption in the unshocked gas to the 
observed levels. Nevertheless weaker flares, generally oc- 
curring at larger radii, may show signatures of cold gas 
absorption as well as "hot" Fe Ka lines. 

Another important aspect of the predicted maximum 
X-ray flare luminosity feauation I38() is that flares occur- 
ring at radii greater than ~ few x W^Rg are too weak to 
be distinguishable on the background X-ray emission of 
Sgr A*. This fact then explains why the flares have du- 
rations in the range from one to tens of kiloseconds (see 

One of the most promising ways to confirm or rule 
out our model is through observations of Fe Ka line emis- 
sion and/or absorption features in X-ray spectra of flares. 
Discovery of such atomic features in Sgr A* spectra would 
provide a very strong support for our model since these 
atomic features are unlikely to be produced in the con- 
text of the three other existing models for Sgr A* flares 
UMarkoff et al. 20011 Liu & Meha 2002; Narayan 2002, 
§5.3). In all of these models flares come from a region very 
close to the black hole, i.e., ^ 10 Rg, because otherwise the 
time scales of the flares would be too long (Baganoff et al. 
2001). The gas (electron's) temperature is very large in 
this inner region, i.e., T ^ 10^^ K in the model by Narayan 
(2002, Fig. 5), T - 10^ - 10" K in Liu & Melia (2002), 
and finally in the jet model of Markoff et al. 2001 the gas 
is relativistic. These temperatures are far too high to al- 
low any noticeable Fe Ka line emission. In addition, gas 
densities are very low in some of these models (e.g., in an 
ADAF or a CDAF), and the recombination time scale for 
completely ionized Fe becomes many orders of magnitude 
longer than the burst duration, again making Fe Ka line 
emission implausible. 

We also found that similar X-ray flares resulting from 
star-disk passages in other LLAGN and nearby inactive 
Galaxies are rather hard to detect except for the most 
nearby and low black hole mass sources. As mentioned in 
fj^l the predicted maximum X-ray flare luminosity (equa- 
tion I38() behaves as Lmax oc l/Mgjj. In contrast, the 
quiescent X-ray luminosity is proportional to oc toMbh 
where m is the dimensionless accretion rate (in the usual 
Eddington units). Hence it becomes more and more dif- 
ficult to discover such flares with increase in the black 
hole mass. Yet some of the nearby sources may show de- 
tectable flares. For example, the recently detected nuclear 
X-ray source in M32 has ~ 0.5 — 5 keV luminosity of about 
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1036 erg/sec (Ho et al. 2003; the black hole mass in M32 
is very similar to that of Sgr A*). This luminosity is low 
enough to allow detection of the strongest X-ray flares, es- 
pecially in the future with instruments such as XEUS and 
Constellation-X. Galaxies not yet detected in X-rays, such 
as M31, hold an even greater promise for future detection 
of X-ray flares. 

Another very promising way to confirm or refute our 
model is through detection of a NIR flare counterpart. In 
our model these flares are emitted some ^ 10"^ — lO'^i?^ off 
Sgr A* location, which translates to ~ O.Of — 0.1 arcsec- 
onds. Current observational techniques do allow to mea- 
sure such a small angular offset to be measured (Genzel 
et al. 2003; Ghez et al. 2003b). Unfortunately the NIR 
emission from flares was predicted to be relatively weak 
(©I although more work is needed to refine this. The 
properties of the disk can also be constrained via stellar 
eclipses and other similar effects (Nayakshin & Sunyaev 
2003; Cuadra et al. 2003). 
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